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i Abstract. In this note we give, for a spectrally negative Levy process, a compact 

formula for the Parisian ruin probability, which is defined by the probability that the 



process exhibits an excursion below zero which length exceeds a certain fixed period 
r. The formula involves only the scale function of the spectrally negative Levy process 
and the distribution of the process at time r. 
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^ . 1 Introduction 
> 
in 

in ' Let X = {Xt,t > 0} be a spectrally negative Levy process on the filtered probability 
^ ■ space {Q,J^,{J^t '■ t > 0},P), i.e. X is a stochastic process issued from the origin which 
^ . has stationary and independent increments and cadlag paths that have no positive jump 

Q ! discontinuities. To avoid degenerate cases, we exlude the case where X has monotone paths. 

As a strong Markov process we shall endow X with probabilities {F^ : x G M} such that 
under Fx we have Xq = x with probability one. Further E^; denotes expectation with respect 
to Pa;. Recall that P = Pq and E = Eq. For background on spectrally negative Levy processes 
^ . we referer to Section 8 of Kyprianou [7j. 



In this paper we deal with the quantity k.^ with r > 0, which is defined by 
Kr = inf{t > r : t — Qt > r} , where gt = sup{0 < s < t : Xg > 0}. 
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Hereby we make the convention that inf = oo and sup = 0. The stopping time is the 
first time the process has stayed below zero for a consecutive period of length greater than r 
and here we are interested in the probability that such an excursion occurs. The number r is 
referred to as the delay. Such Parisian stopping times have been studied by Chesney et al. 
[3] in the context of barrier options in mathematical finance. Dassios and Wu [6] introduced 

in actuarial risk theory: the process X under P^. with a; > is then used as a model for 
the surplus process of an insurance company with initial capital x and the company is said 
to be Parisian ruined if an excursion as described above occurs. We therefore call Kj. the 
Parisian ruin time and the probability of the event {k^ < oo} under P^, the Parisian ruin 
probability. Note that in risk theory, the classical term ruin is referred to as the event that 
the surplus process reaches a strictly negative level. 

In another paper Dassios and Wu [6] gave a formula for the Parisian ruin probability 
and, more generally, the Laplace transform of the distribution of the Parisian ruin time in 
the case where X is a Brownian motion plus drift and in the case where X is the classical 
compound Poisson risk process with exponentially distributed claims. Recently, Czarna and 
Palmowski [1] gave a description of the Parisian ruin probability for a general spectrally 
negative Levy process; hereby they split the analysis into two cases: one in which the 
underlying Levy process has paths of bounded variation and one in which the process has 
paths of unbounded variation. In particular, their result allowed them to reproduce the 
formulas of Dassios and Wu ^ in the two aforementioned cases and to tackle also the case 
where X is the classical compound Poisson risk process perturbed by Brownian motion and 
with exponentially distributed claims. 

Another relevant paper is the one of Landriault et al. [8] . Here the authors study for a 
spectrally negative Levy process of bounded variation a somewhat different type of Parisian 
stopping time in which loosely speaking the deterministic, fixed delay r is replaced by an 
independent exponential random variable with a fixed parameter p > 0. To be a little bit 
more precise, each time the process starts a new excursion below zero, a new independent 
exponential random variable with parameter p is considered and the stopping time of interest, 
let us denote it by Kexp(p) , is defined as the first time when the length of the excursion is bigger 
than the value of the accompanying exponential random variable. Although in insurance 
the stopping time Kexp(p) is arguably less interesting than k^, working with exponentially 
distributed delays allowed the authors to obtain relatively simple expressions for e.g. the 
Laplace transform of Kexp(p) in terms of the so-called (g-)scale functions of X. In order to 
avoid a misunderstanding with the reader, we emphasize that in the definition of Kexp(p) by 
Landriault et al. [8], there is not a single underlying exponential random variable but a 
whole sequence (each one of them attached to a separate excursion below zero); therefore 
IPx('«cxp(p) £ dz) does not equal pe'^'^F^^Hr G dz). 

As the main result of our paper, we provide an expression for the Parisian ruin probality 
which is considerably simpler than the one of Czarna and Palmowski [4j and which simulta- 
neously holds for spectrally negative Levy processes of bounded and unbounded variation. 
Before stating this result, we introduce a little extra notation. 
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The Laplace exponent of X is denoted by 11^(9) , i.e. 



V'(^)=logE[e^^^], 

which is well defined for ^ > 0. We further introduce the scale function W oi X (cf. 
[H Section 8]), which is the strictly increasing, continuous function uniquely defined on 
[0, oo) through its Laplace transform which is given by 

e-'-W{x)<lx = ^ e>0. (1) 

We extend W to the whole real line by setting W{x) = for a; < 0. 

We now give the promised expression for the Parisian ruin probability. Hereby we assume 
that X drifts to infinity (equivalently, E[Xi] > 0) as otherwise Parisian ruin happens with 
probability one. Recall that r stands for a stricty positive real number. 

Theorem 1. Assume E,[Xi] > 0. Then for x > 0, 

r n [fT + z)zF{Xr ^ dz) 

P.K < oo) = 1 - EIX.] J' jL,,p(|,J,^) (2) 

We remark that the methodology of the proof of Theorem [1] can also be used to derive 
the Laplace transform of Kr- However, since the proof would be considerably longer and the 
resulting expression much harder to (e.g. numerically) evaluate, we stick to presenting only 
the Parisian ruin probability. 

In the next section we introduce some additional notation and recall some facts about 
spectrally negative Levy processes which are needed for the proof of Theorem [1] given in Sec- 
tion |3l We conclude by giving some specific examples for which the Parisian ruin probability 
can be expressed in a somewhat more explicit form. 



2 Preliminaries 

We define for a G M, 

r+ = inf{t > : Xi > a}, r" = inf{t > : < a}. 

Note that since is regular for the upper half-line for X (cf. [3 Section 8.1]), we have by 
the strong Markov property r+ = inf{t > : > a}. 
If E[Xi] > (equivalently, ^'(0) > 0), then 

P,(ro-<oo) = l-E[Xi]l^(x), (3) 

cf. [71 Equation (8.7)]. We see in particular that in this case W{x) is bounded by l/E[Xi]. 

The scale function W is absolutely continuous on (0, oo) (cf. [TJ Lemma 8.2]) and we 
denote a version of its density by W. Further, W{0) > if X has paths of bounded variation 
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and W{0) = if X has paths of unbounded variation (cf. [7J Lemma 8.6]). We further recall 
the following well-known expression for the Laplace transform of the first passage time above 
a: 



^<l>(e)(x-a) 



X < a,^ > 0, 



(4) 



where $(^) = sup{A > : ipiX) = 0}, cf. [TJ Section 8.1]. Lastly, we need the following three 
identities. 



Lemma 2. For 9 > 0, 



e-^(0)yW'{x + y)dy, 



m Jo 

POO POO 1 

/ e~'^ / -P(X, G d^)dr = e-*(^)^ y>0, 
Jo Jy r ^{6} 

/ W{z)-V{Xr E dz) =1. 
Jo ^ 

Proof. Using the known identity (cf. [TOj Equation (35)]) 



(5) 
(6) 
(7) 



E., 



' pX - 



e^^' - 7/'(p)eP^' / e-P'W{z)dz 



P 



W{x), p>0, 



by dlD and a change of variables and an integration by parts we get: 



E, 



e ° 1k-<oo} 



.e*W- ( 1 - / e-*(^)^iy(^)d^ 



m 

oo n 

e-*(')^l^(z)dz - ^W{x) 



W{x) 



POO 

--0 / e-*(^)W(x + y)dy 



$(^) 



W{x) 



e 



'^^)yW\x + y)dy. 



m Jo 

This proves the first identity. For the second, we have by Kendall's identity (cf. [H Corollary 
VII.3]), Tonelli and (H, 

/•OO /'OO /*00 /'OO 

/ e-^' / -Po(X, e dz)dr = / e-^' / Po(r+ G dr)d^ 



In order to prove the last indentity, we use again Kendall's indentity, Tonelli and (jl]) combined 
with (II]) to get 

poo POO POO POO 

/ e-^" / W{z)-F{Xr G d2)dr = / e"^" / W{z)¥{t+ G dr)dz 

/•OO 1 

= / e-*(^)W(z)dz = -. 
Jo ^ 



Hence ([7]) follows by Laplace inversion (note that the left hand side of ([7]) is continuous in 
r because a Levy process is continuous in probability and W is bounded and continuous). 
Note that ([7]) can also directly be deduced from [1, Corollary VII. 16]. ■ 



3 Proof of Theorem [T] 

In the case when X has paths of unbounded variation, we will use a limiting argument in 
the proof. For this reason we introduce for e > the stopping time which is defined by 

= inf {t > r : t — gl > r, Xt-r < 0}, where gl = sup{0 < s < t : Xg > e}. 

The stopping time is the first time that an excursion starting when X gets below zero, 
ending before X gets back up to e and of length greater than r, has occurred. Note that 
We have for x < by the strong Markov property and the absence of upward 

jumps, 

< oo) =P,.(r+ > r) + E, [P, (r+ < r, < < oo\J^,+ )] 

=P,(r+ > r) + P,(r+ < r)P,« < oo) (8) 
=l-Px(r+<r)(l-P,«<oo)). 

Using the above we have for x > again by the strong Markov property. 



P,« < oo) =E, 



Px < < ool J" - 



-{to <oo} 



Fx _ « < oo) 



=P,.(ro- < oo) - (1 - Pe« < oo)) E,. 



(9) 



We proceed by finding an expression in terms of the scale function W and the law of Xr 
for the expectation in the right hand side of (Q. By spatial homogeneity, we can write for 

X > 0, 



E. 



lK-<oo|Px_(r,+ <r) 



E, 



[0,oo) 



E, 



[0,oo) 



, ° ^0 



l{r-<oo.-x ^edzj'^oirX, < r) 



Since we have for 6', 2; > by an integration by parts and 



00 1 
e-^^Po(r+ , < r)dr = -Eo 



-$(e)(2+e) 
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we deduce in combination with Tonelli, ([5]) and ([6]), that for 6^ > and x > 0, 

1 



e-^'-E. 







dr 



-$(9)e 
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'S>{e)x _ 



°o p-'i'{^?)(y+^) 



-W(x + t/)d?/ 



POO /^OO POO 

/ W'{x + y) e-^W -Po(X, G dz)drdi/. 

(10) 



Hence by Tonelh and Laplace inversion (noting that both sides of the equation below are 
right-continuous in r), for a; > 0, 



OO PZ — € 



W'{x + y)dy-Fo{Xr G d^) 



J [W{x + z-e)-W{x)]-Fo{Xredz). 



'IV. 



The next step is to prove ([2]) for x = 0. For this we split the analysis into two cases, 
W{0) > and W{0) = 0. First we consider the case 1^(0) > (or equivalently X has paths of 

bounded variation). Then by ([3]), Po(T(7 < oo) < 1 and thus Eq 

Using (E]) and with x = e = and ([3]) and (El), we get: 



< 1. 



Fo{Kr < oo) 



Po(ro- < oo) - /o°"[iy(z) - W{0)]^Fo{Xr E dz) 
I- J,°"[Wiz)-WiO)]^Fo{Xredz) 
^ -W{0)E[X^] + W{0) ^FpjXr e dz) 

Wio) fPo(x. e dz) 

_ f^^Fo{Xredz)-E[X,] 
^Po(X. e dz) 

Now we deal with the second case, i.e. we assume W{0) = 0. Let e > 0. Then F^{tq < 
oo) < 1 (cf. (|3])) and consequently using ([9]) and f lTTj) with x = e and ([3]), we deduce 



PJk! < oo) 



Pe(ro- < oo) - X'"[iy(z) - Wie)]^Fo{Xr E dz) 

1 - nWiz) - Wie)]^n{Xr e dz) 
1/Wie) - E[Xi] + fPo(X. G dz) - r SfPo(X. G dz) 



(12) 



mie) + r 'MXr e dz) - r ^fPo(X. G d.) 
We now want to compute the limit as e J, of both sides of (fT2l) . To this end, recalling that 
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< ]Eo[Xi] < oo, we have by ([7]), an integration by parts and I'Hopital 

40 W{e) 40 Jo ^ 

= 'iff (/' 7^»<^'' ^ *f ) - r ¥^ r ^°<* 

40 W^'(e) 

(13) 

For the hmit as e J, of the left hand side of (|T2l) we introduce for e > the stopping time 

K.^ = mi{t > r : t — gt > r, Xt-r < — e}, where gt = sup{0 < s < t : Xg > 0}. 

We easily see that for < e' < e, {k^ < oo} C {k^' < oo} and Ue>o{K^, < oo} = {kj. < oo}. 
Hence by spatial homogeneity 

limP,« < oo) = limPnffi;!: < oo) = Po(k^ < oo) 
and combined with ( !T2|) and ( !T3|) . this leads to 

Hence, recalling ([7]), we have shown in both cases that ([2]) holds for x = 0. Now plugging 
([3]) and the above expression for F{Kr < oo) into (j9]) and using (|TT1) with e = 0, we arrive at 
(E]) for any a; > 0. 



Remark 3. For x < 0, we have using ([8]) and (|2]), 

P fr^ < 

P.K<oo) = l-E|X,lj^Aji__i_. 

4 Examples 

If one wants to evaluate the formula ([2]), one needs to know the scale function W and the 
distribution of Xj.. There are plenty of examples of spectrally negative Levy processes for 
which the distribution at a fixed time has a closed-form expression and, thanks to recent 
developments (for the latest we refer to Section 3 of Chazal et al. [2]), there are also a lot of 
examples for which the scale function is known explicitly. Unfortunately, it seems that there 
are only a few examples for which both of them are known in closed-form. Below we give 
three examples for which this is the case and for which we give the corresponding expression 
for the Parisian ruin probability. In the other case, one can resort to numerically inverting 
the Laplace transform (which recall is explicity given in terms of ip) of one or both of the 
two quantities. 
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4.1 Brownian motion 

Let Xt = fit+crBt, with fi,a > and {Bt, t > 0} a Brownian motion. Then ip{9) = ii6+\a'^9'^ 
and we easily deduce from ^ that 

Hence by ([2]) for x > 0, 



¥JTr < CX)) =e" 



z¥{Xr e dz) 



Setting X = in above formula, comparing with (fT4|) and realizing that E[Xr] = rE[Xi] = /ir, 
we see that 

/•OO /"OO 

/ e-^^'zF{Xr edz)= / zF{Xr. G dz) - fxr. 
Jo Jo 

Noting that Xr has a normal distribution with mean fir and variance cr^r and making the 
change of variables y = we get 



z¥{Xr e d;^) = / ze ^^dz 

V 27r cr^r Jo 



cr\/r /""^ _y^ fir _^ 

ye dy -\ — = / e dy 



27r J-Hl^ \/2tt J-^tyH. 



OO 5 

2C 



e 2^ + /xrA/" 



72^ V ^ / 

where A/" is the cumulative distribution function of a standard normal random variable. 
Putting everything together results in 

,,2„ 



Fr(Tr < Oo) = C '^t^^JL^ \- X > 0, 



which agrees with the formula found by Dassios and Wu [B]. 

4.2 Classical risk process with exponential claims 

Let Xt = ct- J2i=iCi, where {Nt,t > 0} is a Poisson process with rate t] and Ci, C2, . . . 
are iid exponentially distributed random variables with parameter a independent of A^^. 
Its Laplace exponent is given by ip^O) = c9 — rj + rjg^. Assume that c > rj/a, so that 
E[Xi] = tp'{0) = c-ri/a>0. From it easily follows that 



W(x) = ( 1 - ^e'^-c-^^A , X > 0. 

c — rj/a \ ca / 



Recalling that a sum of iid exponential random variables equals a gamma random variable 
and utilizing the independence between and Nt, we get 



Nr 



J2C^edy \ =5^P^aGdy P(iV, = fc) 

.i=l / A:=0 \i=0 / 

=e--^o(dy) + f:«'=^^^^e--d, 
fe=i ^ ^' 

1 ^-^ ml[m + 1)! 



m=0 



where So{dy) is the Dirac mass at 0. It follows that 

r-oo per / f^,^^\m+l 



znXr e dz) = ze-^^ i^5o{cr - dz) + e''^^-^) ^^^1)! ~ '^"^"^^ 

^ m\{m + 1)! Jo 

m=0 ^ / " u 

oo 



=e I cr 



=e I cr 



m=0 

m\(m + 1)! 

m=0 ^ ' 



crT{m + 1, era) T{m + 2, era) 

a 



where r(a, x) = e ^dt is the incomplete gamma function. Putting everything together 
and using the same trick as in the previous example leads for x > to 

P < oo) ^r-.^. r-nXredz)-ic-v/a)r 



e^^{c — Tj/a) 



c + y^^^" n u 1 1M \cTim + 1, era) — -Tim + 2, era)] 

.!--'m=0 m!(m+l)! L ^ ' ' ar ^ ' 'J 



Although it does not seem easy to show directly that the above expression is equal to the 
one found by Dassios and Wu [6], one can check numerically that the two expressions match. 

4.3 Stable process with index 3/2 

Let Xt = ct + Zti where e > and > 0} is a spectrally negative a-stable process 

with a = 3/2. The Laplace exponent of X is given by ^{B) = c9 + ^3/2. One can m a 
straightforward way check via ([T]) that 

W{x) = ^[l-Ei/2i-c^/^)], x>0, 
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where Eij2{z) is the Mittag-Leffler function of order 1/2, i.e. Ei/2{z) = YlT=o r(ifc+i) • 
follows that for x > 0, 



Prf/tr < OO 



Ei/2{-cy/^T^)zF{Zr edz- cr) 



x\'^r ^ J poo 



zF{Zr edz- cr) 

Using the scahng property of stable processes (see e.g. [U Section Vlll.l]) and the expression 
for the density of the underlying stable distribution taken from Schneider pi] Equation (3.4) 
and (3.5)], we get the following expression for the distribution of Zr'. 



F{Zr G dy) = F{p/^Zi e dy) = { V ^ ' ' ^ ^ ^ 

'^r'/'y-'e-/'W_y2Mu)dy y<0, 



where u = and W^^f^iz) is Whittaker's W-function, see e.g. (9] Section 9.13]. 
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